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Man-made composite materials called "metamaterials" allow for the creation of unusual wave propagation 
behavior. Acoustic and elastic metamaterials in particular, can pave the way for the full control of sound in 
realizing cloaks of invisibility, perfect lenses and much more. In this work we design acousto-elastic surface 
modes that are similar to surface plasmons in metals and on highly conducting surfaces perforated by holes. 
We combine a structure hosting these modes together with a gap material supporting negative modulus and 
collectively producing negative dispersion. By analytical techniques and full-wave simulations we attribute 
the observed behavior to the mass density and bulk modulus being simultaneously negative. 

Classical waves such as sound and light have recently been put to the test in the challenges for cloaking 
objects 1 " 8 and realizing negative refraction 9 " 18 . Those concepts are just a few of recent fascinating phe- 
nomena which are consequences of artificial electromagnetic (EM) or acoustic metamaterial designs. 
Perfect imaging 19 or enhanced transmission of waves in subwavelength apertures 20 are other disciplines within 
the scope of metamaterials which have received considerable attention both from a theoretical and experimental 
point of view 21 " 24 . Smith et al. designed an EM composite material in which the electric and magnetic response act 
simultaneously to exhibit a negative effective index of refraction band 25 . The idea was to tune the resonance of a 
periodic array of split ring resonators such that the negative permeability would occur below the plasma fre- 
quency of a metallic rodded medium where the electric response is negative 26 ' 27 . An advanced version of this 
design is provided by the so-called fishnet structure. In its various build-ups it contains a high figure-of-merit and 
gives rise to negative refraction in the near-IR 10 , microwave 28 and optical regime 13 . 

Manipulating sound waves is readily possible by utilizing metamaterials. Those materials are realized by 
resonating building blocks that are fabricated on a size scale smaller than the wavelength of the irradiated acoustic 
wave. It is the ability to control and tune those meta-atoms, which forms the basics in the design for tailored and 
unusual acoustic material responses. There are different strategies in creating negative refraction using acoustic 
metamaterials. One approach is to use perforations. An isotropic negative index can be obtained by coiling up 
small channels of perforations 17 . Hyperbolic metamaterial for broad-angle negative refraction on the other hand, 
can be constructed using layered holey structures 18 . Another approach is to employ locally resonating structures. 
The first metamaterial fabricated to possess double negativity (simultaneous negative effective bulk modulus 1/k 
and mass density p) was a one- dimensional tube design consisting of periodic interspaced open side branches and 
membranes 29 . Both of those resonators describe a Drude-like behavior for the effective bulk modulus and mass 
density respectively. 

Results 

Designer acousto-elastic surface modes. To go beyond the lumped element design we introduce structured 
metallic components which are able to convey tunable double negativity behavior in a planar configuration 
although extensions into 3 dimensions should not add any difficulties. For the design of such a material we take a 
close look at the EM fishnet metamaterial and initially aim at converting its plasmonic interpretation into 
acousto-elasticity 13 . In its basic design we have two adjacent holey films made out of good conductors (Au, 
Ag) which are separated by a dielectric spacer. Metal surfaces punctured with subwavelength holes allow for the 
EM waves to penetrate much in the same way as a plasmon polariton 30 . Interestingly does the effective 
permittivity take the same plasma form whether they are supported in metals such as silver or in perforated 
highly conducting surfaces. It is this behavior which describes the electric response of a fishnet structure that is 
generated in the holey films. We will start out examining the acoustic analogy of these designer surface modes and 
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Figure 1 | Transmittance and effective mass density spectra of the structured and elastically filled rigid screen. For all simulations, the normalized 
screen thickness is fixed at hld x = 0.7 and the filling material has p = 2300 kg/m 3 and a longitudinal speed of sound q = 1410 m/s. Furthermore we 
consider the sample to be immersed in water, c 0 = 148 1 m/s and p 0 = 1000 kg/m 3 . The transmittance (b) and effective mass density Re(p^) (c) spectra are 
calculated for ajd x = 0.5. In (d,e) we hold c 0 /c t = 2. The full (dotted) lines represent data obtained by modal expansions (COMSOL simulation). The 
incident sound plane wave is impinging at the normal direction, (f) The band diagram is simulated for a structure with the parameters ajd x = 0.5 and hld x 
= 0.7 and contains filling parameters as in the latter example though with c 0 /c t =3.7. 



later study the effective bulk modulus which is showing resemblance 
to the magnetic response of the fishnet structure. 

Let us suppose that we have a perfect rigid screen into which holes 
are carved. In order for an incident airborne sound field to be highly 



localized, we fill these indentations with an elastic material as ren- 
dered in Fig. la. We begin with a simple model that is constructed for 
a structure with translational invariance in the y direction and tex- 
tured by periodic slits of lattice constant d x along the x axis. An 
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incident wave excites elastic cavity modes inside the filled indenta- 
tions which predominantly are governed by a displacement u z in the 
direction perpendicular to the lattice. The displacement generally 
satisfies the homogeneous wave equation for an isotropic solid: 



d 2 u z 



fd 2 u z 



d 2 u z \ /n Jd 2 u z 



dxdz 



(i) 



where k, \i and p are the modulus of incompressibility (first Lame 
coefficient), modulus of rigidity (second Lame coefficient) and the 
solid mass density respectively. As the elastic inclusions are clamped 
to the rigid frame, u z (x = 0, a x ) = 0 must vanish at the edges of the 
slit. Furthermore, we must guarantee zero shear- stress at the free- 
space and metamaterial interface as waves only propagate longit- 
udinally in a fluid, o xz = 0. Given the long wavelength limit which 
we are mainly interested in, i.e., k > a x , we take the fundamental 
cavity eigenmode to be dominant as it penetrates the deepest into the 
elastic inclusion, hence, we can write the out-of-plane displacement 
field inside this region as: 

u z (x,z) =A sm(P z z)' sin (^-x^j -e iwt . (2) 

Disregarding in -plane displacement (u x = 0) is a plausible assump- 
tion in narrow elastic channels, hence we substitute Eq. (2) into Eq. 
(1), whereupon we solve for the wavenumber /? z comprising longit- 
udinal and transverse wave motion: 



P z = co 



k + 2fi 



p a x co 1 c\ 



(3) 



which possesses a cut-off frequency co p = c t — .As we are focusing on 

the long wavelength regime we take into consideration only specular 
reflection when imposing continuity of an normally incident plane 
wave irradiating onto the structured half-space (h — > °°). From here it 
is straightforward to write down expressions for the reflection co- 
efficient and the effective impedance: 



8 a x 



yC 



(4) 



where the structure factor Sf = —z-r- Substituting Eq. (3) into Eq. (4) 

n l d x 

and rewriting the effective speed of sound as c e ff = k 0 /p z , we arrive at: 



Peff=- 



K 2 d x pco 



when solving for p e p which is the canonical Drude form of an 
acousto -elastic surface mode with a geometrical "plasma" frequency 
of 



co p = c t - 



(6) 



the cut-off frequency of the elastic waveguide. In summary we can say 
that structured rigid screens will have this form of the effective mass 
density provided all dimensions are treated in the effective medium 
limit. The acousto -elastic cut-off frequency can be modified at will 
such that the screen can constitute a tunable high-pass filter with a 
stopband from zero frequency up the cut-off. 

In the following we will take a look at the tunability of the struc- 
tured screen with the aim to design the spectral location of the 
acoustic "plasma" frequency. In order to do so, we employ a rigid 
structured screen of finite thickness h as depicted in Fig. la and 
calculate complex transmission and reflection coefficients for using 
them to numerically obtain the effective constitutive parameters 



based on a retrieval technique 31 . The theoretical formalism used to 
calculate the scattering coefficients is based on the modal expansion 
of the pressure and velocity fields in free- space such as the displace- 
ment and stresses inside the elastic inclusion 32 . Even though, in the 
former, the metamaterial halfspace is treated in the effective medium 
limit (k >> a x ) by solely introducing specular reflections and taking 
into account only the fundamental elastic cavity mode, we conduct 
simulations by considering all higher expansion orders. As given in 
the caption of Fig. 1, we hold the screen thickness h/d x that is normal- 
ized to the lattice d x fixed, but vary the width a x ld x of the indentations 
and also the inclusion material by different values of the transversal 
speed of sound c t . Transmittance spectra are plotted and compared to 
the real part of the effective mass density Re(p e ff) as depicted in 
Fig. lb and Fig. lc respectively. We clearly see that by lowering c t 
we also lower the onset of the first allowed band which means that 
the spectral locations of the acousto -elastic plasma frequencies 
shift towards longer wavelengths kldx. These are the locations at 
_ ^Jtco _ ^ Co yfaexe the effective mass density, as seen in 

(D p C t 

Fig. lc, changes sign - the transition between an opaque and trans- 
parent effective fluid. Fig. Id and Fig. le confirm that the acousto - 
elastic plasma oscillation can be further controlled by the indentation 
width a x . With a specific chosen filling material, in this case c 0 /c t = 
2.0, we predict that the cut-off wavelengths now locate at k p = 4a x . In 
order to verify the analytical approach based on modal expansions 
(full lines in Fig. 1), we have further conducted full- wave simulations 
(COMSOL Multiphysics, dotted lines) and found very good agree- 
ment. In addition we calculate the radiative transmittance band dia- 
gram which has been plotted in units of the plasma frequency. This is 
rendered in Fig. If (see the captions for details) where the low fre- 
quency plasmon-like mode interestingly shows very little dispersion 
with the parallel momentum k x of the incident sound. Again it must 
be noted that the cut-off can be controlled at will by the size of the 
indentations and the choice of the filling material such as polymers. 
This allows for the creation of designer surface states at almost arbit- 
rary frequency. 

Tunable negative effective bulk modulus. The magnetic response of 
the aforesaid EM fishnet structure is induced by the component 
consisting of metal-insulator-metal (MIM) units. These structures 
are hosting resonances which can give rise to negative values of the 
effective permeability. Recently a theoretical study on the acoustical 
version of the MIM structure showed that these sandwich-materials 
would give rise to an effective negative bulk modulus 1/tc, which su- 
stains for a broad range of angles of incident sound 33 . The conse- 
quence of a negative bulk modulus in a fluid element is the overall 
expansion of it as a reaction to a positive acoustic stimulus. This 
behavior has recently been observed experimentally in the audio-fre- 
quency range 34 . In Fig. 2a we show the slitted version of this design 
and compute the normally incident sound transmittance through it 
by the same method as reported in 33 . We fill a fluid in the interme- 
diate gap region of the structure (with height h g ) and study the con- 
trollability of the resonance by changing the value of the speed of 
sound Cg. We compare the spectral transmittance (Fig. 2b) to the 
retrieved effective bulk modulus Re(l//c) (Fig. 2c) that at resonance 
sustain negative values over extended ranges. These regions of nega- 
tive 1/k are as expected accompanied by polaritoniclike bandgaps, 
which in Fig. 2b are seen by the presence of sharp dips in the 
transmittance spectra. Moreover, we predict a distinct correlation 
between these gap induced resonances and the speed of sound c g 
therein as the resonating modes shift towards lower frequencies 
when the filling material speed-of-sound is lowered. 

Composite material exhibiting double negativity. The aforemen- 
tioned tuning strategy will prove to be very useful at the stage when 
merging the previous Drude-like metamaterial containing negative 
/^together with the present one to form a composite demonstrating 
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Figure 2 | Transmittance and effective bulk modulus Re(l/A*) spectra for 
two facing structured rigid plates with a/d x = 0.5, hld x = 0.2 and the gap 
separation h g ld x =0.1. The gap is filled with a fluid of relative mass density 
Pg/po = 1 and speeds of sound cjc 0 as indicated. The transmittance is 
presented in logarithmic scale and calculated at normal incidence. 



double negativity. The pursuit is to shift these bulk modulus 1/k 
resonances shown in Fig. 2, below the cut-off frequency co p of the 
elastomer filled screen. This particular design idea is a utile extension 
of the former one, we solely need to fill the slits with an elastic 
material as sketched in Fig. 3a. We design a plasma frequency to 
be located at k p ~ 2.%d x below which the effective mass density is 
negative. We have chosen to use a fluid inside the gap of speed c g /c 0 = 
0.25 such that the negative response of 1/k occurs within the 
forbidden region, below the plasma frequency. The combination of 
simultaneous negative constitutive parameters induces a passband of 
full sound transmission at around X ~ 4.4<i x , as seen in Fig. 3b. Here 
again we compare analytical predictions with full-wave simulations 
where a fairly good agreement is obtained, confirming propagation at 
a region of double negativity. Interestingly we observe negative 
dispersion in this propagation band which has been computed by 
the transmittance as a function of frequency and parallel momen- 
tum, as rendered in Fig. 3c. 

Discussion 

We have demonstrated both by analytical and numerical simulations 
that structured rigid screens filled with elastic inclusions can be 
described in the long wavelength limit as an effective medium char- 
acterized by a mass density of the plasmon form. This type of meta- 
material is thus not only used to control surface waves but also to 
create them. The combination of this behavior together with gap 
materials hosting bulk modulus resonances, produces negative disper- 
sion with simultaneous negative values of p e ff and l//c ej ^ which is a 
characteristic of negative refractive index materials. The near-field 
coupling of the constituting resonances can be considered weak, which 
is essential when tuning the two effective parameters independently. 



0.50 




-10 



3.8 



4.0 



4.2 



4.4 

X/d 



4.6 



0.25 



V 



(b) 


41 i 

h ;; h g /d x" 01 

1 : • h/d = 0.5 
I* : 

If •• c /c = 0.25 
Is • ? 0 




w 

• 1 

' 1 \ 



1.0 



0.8 



0-6 g 



0.4 § 
3 



0.2 



0.0 



5.0 



h/d = 0.1 

8 x 

h/d x = 0.5 
c/c =0.25 



0.25 

kd/Ti 



0.50 



Figure 3 | Acoustic double negativity band in a composite metamaterial. (a) The composite metamaterial is designed with parameters ald x = 0.5, hld x = 
0.5 and hjd x =0.1 with an elastic inclusion containing p = 2300 kg/m 3 , c\ = 1400 m/s and c t = 530 m/s. (b) With Cg/c 0 = 0.25 the transmittance spectra 
(right axis) is accompanied with plots of Re(l/K: eJ ^) and Re(p e ^) (left axis). The full (dotted) lines represents data obtained by modal expansions 
(COMSOL simulation), (c) Band diagram simulation for the structure with parameters as given in (b). 
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Methods 

We have used modal expansion techniques (rigorously coupled modes theory) for 
analytically deriving the reflection coefficient and the impedance leading to the 
expression of the effective mass density in perforated elastically filled screens, Eq. (5). 
In the second stage, this technique has been implemented to numerically compute the 
transmission data both as a function of frequency co and parallel momentum k x 32 ' 33 . 

Based on an effective medium retrieval technique, we were able to calculate the 
effective acoustic paramaters such as the mass density p ej ^and the bulk modulus II 
We have employed the method designed by D.R. Smith et al. and slightly modified it 
for the case of sound 25 . 

The finite element method program COMSOL Multiphysics was utilized to con- 
duct verifying numerical experiments seen in Fig. lb - Fig. le, such as Fig. 3b. 
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